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Abstract 

In  this  effort  we  consider  damage  detection  algorithms  using  thermal  probes. 
The  techniques  proposed  employ  models  derived  from  homogenization  theory  as 
the  model  solution  in  parameter  estimation  procedures  for  simulated  data  for  heat 
flow  in  a  porous  medium.  We  consider  data  simulated  from  a  model  on  a  perforated 
domain  with  isotropic  flow  and  data  simulated  from  a  model  on  a  homogeneous  do¬ 
main  with  anisotropic  flow.  We  report  on  ordinary  least  squares  methods  combined 
with  the  method  of  mappings  to  detect  the  presence  of  damage  and  to  characterize 
the  geometry  of  the  damage. 
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1  Introduction 


Nondestructive  evaluation  techniques  play  an  important  role  in  safely  using  many  struc¬ 
tures  as  they  age  beyond  their  design  lives;  these  structures  include  bridges,  pipes,  and  air 
craft.  Though  nondestructive  evaluation  has  been  studied  and  used  for  decades,  until  re¬ 
cently  much  of  the  methodological  development  has  been  for  homogeneous  materials.  As 
composite  materials  have  become  more  popular,  there  is  a  need  to  develop  nondestructive 
evaluation  techniques  that  accommodate  the  heterogeneous  structure  of  these  materials. 
A  common  thermal  nondestructive  evaluation  technique  is  the  flash  heat  experiment,  the 
use  and  underlying  theory  of  which  has  been  thoroughly  developed  for  homogeneous  ma¬ 
terials  [11,  15,  25].  The  largest  order  heterogeneity  in  composite  materials  is  their  porosity 
which  is  associated  with  the  manufacturing  process.  In  this  effort,  we  include  this  poros¬ 
ity  in  simulated  thermal  data  and  subsequent  analysis  of  our  models.  We  build  on  our 
previous  work  [3]— [5]  which  relies  heavily  on  results  in  [6,  7,  9,  14,  16,  18,  20,  21,  23,  26]. 
In  Section  2  we  discuss  two  models  that  we  use  to  simulate  data;  we  then  consider  damage 
detection  in  Section  3  in  the  context  of  statistical  model  comparison  problems,  and  finally 
we  use  ordinary  least  squares  parameter  estimation  to  consider  damage  characterization 
in  Section  4. 


2  Mathematical  and  statistical  models 

In  order  to  present  our  method  of  detecting  and  characterizing  damage,  we  discuss  a  few 
mathematical  and  statistical  models.  We  developed  a  method  for  modeling  the  flash  heat 
experiment  on  both  damaged  and  undamaged  porous  domains  in  [3]  which  we  use  here 
to  simulate  data.  We  model  the  undamaged  domain  as  a  randomly  perforated  domain  fl. 
The  corresponding  undamaged,  homogeneous,  non  perforated  domain  is  denoted  by  f!  and 
here  is  taken  as  an  L\  x  L2  rectangle  O  =  [0,  Li]  x  [0,  L2]  (we  will  use  L\  =  10  and  L2  =  2). 
We  assume  nr  randomly  placed  pores,  generated  using  methods  described  in  [3],  which 
are  denoted  flj  with  boundaries  dVti  for  i  —  1 ,  2  . . .  nr .  The  undamaged  perforated  domain 
fl  is  given  by  (1  \  (U]!l1flj).  The  four  boundaries  of  fl  which  are  also  the  four  exterior 
boundaries  of  f!  are  denoted  oj,  for  i  =  1,2,  3, 4  as  depicted  in  Figure  1.  We  model  the 
flash  heat  experiment  which  approximates  an  experiment  where  the  bottom  boundary 
uq  =  {(x,y)\x  G  [0,Li\,y  =  0}  is  heated  by  a  flash  heat  source  [25].  Throughout  this 
document,  we  will  refer  to  uq  as  the  source  boundary. 

We  denote  by  £lD(q)  the  damaged  porous  domain  which  is  constructed  using  the  ellipse 

with  center  (x_d,2/d),  horizontal  semi-axis  ho  and  vertical  semi-axis  Vd •  The  damage 
(assumed  to  be  in  the  shape  of  an  elliptical  section)  and  the  damaged  domain  are  pa¬ 
rameterized  by  q  which  specifies  (xD,yD,hD,vr>)  although  throughout  this  presentation 
we  simplify  (without  loss  of  generality)  our  investigations  by  considering  y^  =  2.25  to 
be  fixed.  We  restrict  q  so  that  the  center  of  ED(q)  is  above  the  back  boundary  oj2  of  Cl, 
and  ED(q )  intersects  the  boundary  of  Cl  at  two  points  on  uq.  The  damaged  geometry  is 


2 


(jU2 


Figure  1:  An  example  randomly  perforated  domain  (enlarged  view). 


then  given  by  flD(q)  =  (O  U  O(q))  \  ED(q)  where  O(q)  is  the  set  of  ellipses  hh  in  VI  that 
intersect  ED(q).  We  denote  (with  possible  reordering)  the  boundaries  of  the  n®(q)  pores 
remaining  in  flD(q)  by  dQiy  i  =  1, 2, . . .  n®(q).  The  back  boundary,  c u-f(g),  of  &D(q)  is 
parameterized  by  (x_d,  Ud,  hDl  ud).  In  order  to  define  cyf(g)  more  precisely,  we  introduce 
the  piecewise  defined  function 


To 


r(x]  q)  =  < 


Vd  -  vD\jl 
L2 


X  —  X£> 


h 


D 


for  0  <  x  <  x^q) 

2 

for  XL^q)  <  x  <  xR(q) 
for  xR(q )  <x<Li, 


(1) 


where  the  left  end  point  x^q)  is  given  by  x^q)  =  Xr  —  hn\  1  — 


Vd  ~  L-2 

vd 


and  the 


right  end  point  is  given  by  xR(q)  =  xd  +  hp\  1  — 


Vd  ~  T; 
vd 


The  back  boundary  of 


ffD(g)  is  then  given  by  cuf  (g)  =  {(a;,  y)  :  0  <  x  <  Lx  and  y  =  r(x ;  q)}. 

An  example  damaged  domain  is  depicted  in  Figure  2.  We  model  the  heat  equation  on 
this  domain  with 

A";"1  -  av  •  (VtiJT1)  =0  in  JiD(9)  x  (0,  T) 


9td 

_  „  rand  n 

O'  ^  ^  —  U 


dy 

d 


_  ”  „  rand  n 
Q?  U  rj  —  U 


difD 
a^rjU'“d  = 

ura-nd(0,x)  =  0, 


on  x  (0,  T ) 

on  (cji  Uwf  (q)  U  cu3)  x  (0,  T) 
on  c j4  x  (0,  T) 


(2) 


with  thermal  diffusivity  a  =  2.9167,  flash  heat  intensity  Sf  =  2.75,  and  flash  heat  duration 
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Figure  2:  An  example  geometry  QD(q)  for  (xd,2/d,^c >,Vd)  =  (5,2.25,2,1) 


ts  =  0.6.  The  heaviside  function  is  given  by 


%*.](*) 


1,  for  t  G  [0,  ts] 
0,  otherwise 


We  would  like  to  model  data  that  would  result  from  the  output  of  an  IR  camera.  Based 
on  experimental  protocols,  we  assume  that  each  pixel  gives  the  average  temperature  across 
a  length  i  on  the  source  boundary  0J4.  We  denote  the  corresponding  noise  free  observations 

by  ^ 

(u^nd{q))..  =  jf  s,0-,q)ds,  (3) 

J  Xi 

with  pixel  length  £  =  0.57.  To  include  error  associated  with  the  measurement  process,  we 
will  take  data  given  by  realizations  of 

D“"d  =  («'D“d(4o))«  +  0o/%,  (4) 

where  q0  is  the  usual  “true”  parameter  value  in  statistical  models  [6] ,  the  /3tJ  are  indepen¬ 
dently  identically  distributed  standard  normal  random  variables.  That  is  /3tJ  A/"(0,1) 
(we  use  boldface  to  denote  random  variables  throughout).  We  will  use  00  =  0.015  in 
all  our  calculations.  It  is  important  to  emphasize  that  D^nd  is  a  random  variable  with 
realizations  D™d.  The  realization  D™nd  would  correspond  to  an  observed  data  set.  We 
take  m  spatial  nodes  which  correspond  to  each  pixel  that  begins  at  {ay}  in  (3)  and  n 
temporal  nodes  which  correspond  to  the  in  (3). 

Based  on  the  encouraging  results  in  [5]  which  verified  that  a  homogenization  approach 
provides  good  approximations  to  solutions  urand  for  problems  on  porous  domains,  we  will 
use  a  model  motivated  by  homogenization  theory  as  the  model  solution  in  our  inverse 
problems.  We  will  also  simulate  data  from  this  model  as  in  [5]  to  understand  the  effect  of 
the  error  (as  compared  to  data  simulated  with  w^nd)  associated  with  the  approximation 
derived  from  homogenization  theory  on  the  inverse  problem.  Methods  used  in  [2,  4,  5,  14, 
16,  17,  18,  19,  20,  21,  24]  can  be  used  to  establish  the  good  approximation  of  uT^nd(t,x;  q) 
by  the  homogenization  solution  uD{t,  x;  q)  where  uD(t,  x\  q)  is  given  by 
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(5) 


'pvv^JD  -  «V  •  ( A°VuD )  =  0 

in  r  X  ClD(q) 

a— - uD  =  0 

OljAO 

on  T  x  (a Ji  U  c l>2  (q)  U  u>3) 

on  T  xu4 

uD(  0,x)  =  0 

for  x  G  £lD, 

with  thermal  diffusivity  a,  proportion  of  the  material  in  the  porous  domain  py  (here, 
pv  =  0.98),  anisotropy  matrix  associated  with  the  results  of  homogenization  theory  A0  (see 
[4]  and  [5]),  intensity  of  the  heat  lamp  Sf  and  heaviside  function  Zjo ,ts\(t)-  The  geometry 
on  which  the  partial  differential  equation  (5)  is  dehned  is  given  by  ClD(q)  =  Cl  \  ED(q). 

The  weak  solution  of  (5)  is  given  by 


Pv 


d_ 

dt 


uD(t,  cc ;  q)4>(x)  dx 


+ 


a  /  \/uD(t,x;q)  ■  A°W(j)(x)  dx 

Jn°(q) 


(0)  ^0)4  1 


(6) 


for  all  test  functions  0  G  ih1(fiD)  where  yW4  is  the  trace  on  u4  and  £W4  is  a  parametrization 
of  uj4.  In  order  to  compute  a  finite  element  solution  of  (6),  for  each  iteration  of  q ,  a 
new  parameter  dependent  mesh  of  ClD(q)  must  be  generated.  This  presents  significant 
computational  challenges  in  carrying  out  inverse  problems.  In  order  to  avoid  this  difficulty, 
we  employ  ideas  from  [10,  11]  and  transform  the  homogenized  damaged  domain  ClD(q) 
to  a  parameter  independent  reference  domain  (in  this  case  fl)  where  we  can  readily  use 
parameter  independent  finite  elements  to  carry  out  the  computations.  The  coordinate 
transformation  is  embodied  in  x  =  T(q)  o  z  where  T(q)  is  a  bijection  that  maps  to 
ClD(q)  given  by 


x  = 


x 


—  T(q)(z)  — 


z  1 

rjzy  q ) 
Lo 


-Z2 


(7) 


The  weak  formulation  (6)  can  then  be  written  as  (here 


is  the  transpose) 


Pv 


u(t,  z ;  q)4>(z)  det  (VT(z;  q))  dz 


(8) 


-a  Vu(f,  z\  q)  ■  VT(z-  q)~lA 0  (VT(f;  g)f)  1  V(j>(z)  det  (VT(z;  q))  dz 


in 


=  Sfl[0,ts](t)  /  7W4 (</>)  d&j 4,  for  all  0  G  H1^), 


'  U)4 
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by  making  the  coordinate  transformation  x  —  T(z]  q )  with  u  on  12  dehned  by  u(t,  z\  q)  = 
uD(t,T(z;q);q).  Noting  that 


VT(z-q)  = 


1  0 

r'(zi,q)  r(zi,q ) 
Lo  Z 2  Lo 


(9) 


and  det(VT(£;  q))  =  -  y-  ^  ,  we  may  rewrite  (8)  as 

L2 

T~  f  U  g)0(^)r-(zi;  q)dz  -  a  f  Vu(t,  z;  q)  ■  B(z ;  qr)V0(z) 

^2  Jo  ot  Jn 


'  UJ4 


=  S/Z[0it<](t)  /  7^4(0)  ^4- 
In  (10),  5(z;  g)  is  a  2  x  2  matrix  with  entries 
n  r0b  9) 


fen  (2;  9)  =  an-  T 

Xi  2 

fei2(X;  9)  =  -a^  z2 


fe2i(X;  9)  =  o-q  z2 


0  „  9)  ,  0 


l2 

0  „  9)  ,  0 


u 


+  o12 

+  a21 


h  (?.n)  _  7o  ,  v2(r'(zl]q))2  r\Zl-q ) 

WM,-“,iW  r(„;,)L2  -7(77) 


^2(«12  +  02i)  +  0° 


L2 

22r(z1]qy 
~.Nr 


(10) 


(11) 


where  a^-  are  the  entries  of  Au.  Using  the  finite  element  approximation  uiy(t,z]q)  = 

N 

uf  (t;  q)(j)f  jz)  where  the  (f>f  (z)  are  piecewise  affine  two  dimensional  basis  elements 
3= 1 

(e.g.,  see  [12]),  we  may  solve  for  the  time  dependent  coefficients  uN  =  (u^, . . . ,  w^)l  by 
solving  the  ordinary  differential  equation 


where 


PvCN(q)jfUN{t )  +  aI<N{q)uN(t )  =  S fl[0!ts](t)  f* , 


W  =  L 


2  Jn 


<Pj  {z)<t>i  ( z)r(Zl;q)dz  for  ij  =  1,2, ...  N 


(12) 


KS{q)  = 


( z )  •  B{z;  q)  V0f  (z)  dz  for  i,  j  =  1, 2, ...  N 


(13) 


fi 


PN 


iuMi  )di 


for  i  —  1,2 , ...  IV. 


'  CJ4 
rand/ 


As  in  formulating  (ur£n  (q))ij  in  (3),  we  use  the  average  of  the  solution  over  intervals 
of  length  £  =  0.57  (which  represents  the  pixel  length)  in  our  model  solution.  We  will  use 
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model  solutions  given  by  (recall  z\  =  x) 

Uiji(l)  =  \  f  uN(tj,x,0;q)dx  (14) 

where  uN(t,  z;  q )  is  the  finite  element  solution  obtained  using  (12). 

We  will  use  UjJ(q)  as  defined  in  (14)  as  a  model  solution  in  our  data  analysis.  To 
consider  the  effect  of  model  discrepancy  on  our  statistical  procedures,  as  in  [5],  we  also 
consider  (in  addition  to  the  data  generated  using  (4))  data  generated  using  the  finite 
element  solution  given  through  (12).  Specifically,  we  will  consider  data  generated  using 

D«  =  Dg  =  E#(®,)  + ,70ft,.,  (15) 

where  as  in  (4),  f3t]  ~  J\f( 0, 1).  The  realizations  of  the  random  variables  D ^  are  denoted 
D ij.  We  remark  that  we  will  suppress  the  index  N  since  it  will  be  fixed  throughout. 
Computational  studies  of  choosing  appropriate  values  of  N  in  (12)  and  of  ND  in  corre¬ 
sponding  finite  element  solutions  of  uT^nd  in  (2)  were  carried  out  before  beginning  our 
inverse  problem  calculations.  Values  for  N  ranging  between  48  and  2337  and  for  ND 
ranging  between  64  and  3193  were  tested.  Above  N  =  609  and  ND  =  829  there  was  little 
change  in  computed  solutions  although  computational  times  increased  substantially  for 
larger  values  than  these.  Therefore  we  chose  these  values  for  the  computational  results 
for  detection  and  characterization  reported  on  here. 

We  will  consider  two  types  of  statistically-based  procedures.  We  will  formulate  a 
comparison  technique  to  detect  damage  using  a  model  comparison  test  which  is  devel¬ 
oped  in  [8]  and  illustrated  in  [6,  12].  We  will  then  demonstrate  use  of  ordinary  least 
squares  parameter  estimation  to  estimate  the  parameters  (xo,  ho,  Vo)  which  characterize 
the  damage. 

3  Damage  Detection 

Both  the  damage  detection  methodology  and  the  ordinary  least  squares  (OLS)  parameter 
estimation  involve  minimizing  the  OLS  cost  functional  which  is  given  by 

m  n 

J(9)  =  EE  (Uij(q)  ~  d&t&ij)2.  (16) 

i= 1  j= 1 

(As  noted  above,  since  N  is  assumed  fixed  throughout,  we  will  suppress  the  corresponding 
superscripts  on  variables  such  as  U*  (q)  and  numerous  others  throughout  the  remainder 
of  this  presentation.)  The  term  data^  represents  without  loss  of  generality  either  of  the 
realizations  or  D)“ld.  Here  we  emphasize  that  the  function  J(q)  of  the  realization 
data*j  is  a  realization  of  a  function  of  a  random  variable  datajj  given  by 

m  n 

J(v>  =  (UM  ~ datau)2-  (17) 

i=i  j= i 


7 


where  data*,  represents  either  the  finite  element  based  random  process  or  the  random 
porous  domain  random  process  DDnd. 

Damage  detection  may  be  considered  using  a  hypothesis  or  model  comparison  test. 
The  case  where  there  is  no  damage  can  be  modeled  with  an  ellipse  ED(q )  that  has  a 
horizontal  semi-axis  length  of  zero.  So  when  hp  =  0,  we  have  that  there  is  no  damage 
in  the  sample.  We  use  this  to  formulate  a  hypothesis  test.  The  null  hypothesis  (there 
is  no  damage)  is  given  by  H0  :  ho  =  0,  and  the  alternative  hypothesis  (there  is  possible 
damage)  is  given  by  Ha  :  ho  %  0.  We  will  define  the  set  of  parameters  corresponding  to 
the  null  hypothesis  by  Qh  —  {q  —  {x^,  ho,  Ud,  ot)\hr>  =  0  ,q  E  Q}.  The  full  admissible 
parameter  set  Q  —  {q  —  (xD,hD,V£,,a)  e  R4}  is  assumed  to  be  a  bounded  subset  of 
the  nonnegative  quadrant  of  R4.  Using  methods  described  in  [6,  8,  12],  we  consider  two 
parameter  estimates 

qn  =  arg  min  J(q),  (18) 

q&O-H 


and 


q  =  arg  min  J(q), 
q£Q 


(19) 


where  the  parameter  estimates  are  realizations  of  random  variables.  The  corresponding 
two  estimators  are  given  by 

qn  =  arg  min  J(g),  (20) 

q&Qn 


and 


q  =  argmin  J(g), 
q&Q 


(21) 


respectively.  The  estimates  are,  of  course,  realizations  of  the  random  estimators. 
This  model  comparison  methodology  uses  the  test  statistic  X  given  by 


X  =  nm 


(  J(qn) 

\  m  ) 


(22) 


where  nm  is  the  total  number  of  observations  (there  are  n  temporal  nodes  and  m  spatial 
nodes).  The  random  variable  X  then  has  realizations  X  which  are  given  by 


nm 


(  J(qh )  -  J(q)\ 

\  m  )  ’ 


(23) 


where  qn  is  the  corresponding  realization  of  qn  and  q  is  the  corresponding  realization  of 

Q 

The  results  of  [8]  indicate  that  if  H0  is  true  then  the  random  variable  X  approaches  in 
distribution  (as  sample  size  nm  increases)  a  random  variable  which  follows  a  chi-square 
distribution  with  one  degree  of  freedom  (denoted  y2(l)).  This  can  be  used  in  the  following 
way.  If  we  would  like  to  reject  our  null  hypothesis  with  95%  confidence,  we  would  calculate 
the  corresponding  significance  level  a  =  0.05  (note  that  confidence  =  (1  —  a)  x  100%)  and 
the  corresponding  threshold  value  r  defined  by  the  probability  P(X  >  r)  =  a.  Using  a 
y2(l)  cumulative  distribution  function  table,  for  a  =  0.05,  we  find  r  =  3.84.  Then  for  any 
realization  X  with  X  >  3.84,  we  may  reject  with  at  least  95%  confidence.  The  minimum 


Pixel  Set 

Xi  values 

Left  11  Nodes 

Xi  =  £{i  —  1)  for  i  e  {1,  2, . . .  11} 

Right  11  Nodes 

Xi  =  £(i  —  1)  +  6 £  for  i  E  {1,2,...  11} 

Evenly  Spaced  11  Nodes 

Xi  e  {0, 2£,  5£,  6£,  7£,  9£,  10£,  12£,  13£,  15£,  16£} 

All  17  Nodes 

Xi  =  £(i  —  1),  for  i  E  {1,2,3...,  17} 

Table  1:  Pixel  sets  used  throughout  the  computations  with  £  =  0.57 


value  a*  at  which  the  null  hypothesis  can  be  rejected  for  a  realization  X  is  called  the 
p- value,  that  is 

p- value  =  a*  =  min P(X  >  A"). 

Thus,  the  closer  a  p-value  is  to  zero,  the  more  confidence  with  which  one  may  reject  the 
null  hypothesis  (for  further  discussion  see  [6,  8,  12]). 

We  considered  simulations  using  four  different  sets  of  spatial  nodes  (corresponding  to 
pixels)  which  are  listed  in  Table  1.  We  used  these  values  to  generate  the  data  Dy  and 
DUnd.  We  considered  eighteen  different  damages.  In  each  example  we  assumed  yo  =  2.25 
is  known.  We  considered  one  series  of  nine  damages  of  increasing  size  with  center  xd  =  3 
and  another  series  of  nine  damages  increasing  size  with  center  xd  =  5.  For  these  examples, 
we  used  vlD  =  0.3  +  0.05(i  —  1)  and  hlD  =  2vlD  for  j  G  {1,2,3... 9}.  Example  damages  for 
xjj  =  3  and  i  =  1,2,8,  and  9  are  depicted  in  Figures  3  (a),  (b),  (c)  and  (d),  respectively. 
In  Figures  4-7,  we  depict  for  comparison  the  results  of  this  damage  detection  technique  for 

data  sets  simulated  with  (4)  and  (15)  using  spatial  nodes  given  in  Table  1  and  temporal 

1 

nodes  given  by  tj  =  0.6  +  f°r  3  £  {1,2,...  92}. 

The  threshold  value  of  the  statistic  r  at  which  one  may  reject  the  null  hypothesis 
with  99.9%  confidence  is  r  =  10.8.  Thus,  if  X  >  10.8  one  may  reject  the  null  hypothesis 
with  99.9%  confidence.  The  larger  X,  the  more  confidence  there  is  in  rejecting  the  null 
hypothesis. 

In  Figures  4-7  (a)  we  depict  the  test  statistic  X  using  (23)  for  xd  =  3,  and  in  Figures  4- 
7  (b)  we  depict  the  test  statistic  X  using  (23)  for  xd  =  5,  in  Figures  4-7  (c)  we  depict  the 
p-value  for  xd  =  3,  and  in  Figures  4-7  (d)  we  depict  the  p-value  for  xd  =  5.  We  observe 
that  for  the  damages  above  3%  the  test  statistics  X  become  very  large.  There  is  no 
clear  difference  between  the  values  of  the  test  statistic  X  when  calculated  using  D ^  than 
when  using  D™nd  which  suggests  that  the  model  discrepancy  associated  with  the  difference 
between  ur£nd  (the  solution  of  (2))  and  uD,  the  solution  of  (5),  or  its  approximation  given 
by  (12),  does  not  substantially  affect  this  test  statistic. 

Using  the  left  most  nodes  (see  Table  1)  in  Figures  4  (c)-(d),  Figures  5  (c)  and  (d), 
we  depict  the  p- values  for  xd  =  3  and  xd  =  5,  respectively.  For  both  of  these  series 
of  damages,  the  p-values  calculated  using  Dj“ld  are  very  low  suggesting  that  one  may 
reject  the  null  hypothesis  with  high  confidence.  The  p-values  calculated  using  D,,,-  are 
larger,  though  still  sufficiently  small  to  reject  the  null  hypothesis  with  90%  confidence. 
In  every  example,  for  damages  which  occupied  more  than  3%  of  the  area,  the  p-values 
were  close  to  zero  which  corresponds  to  high  confidence  in  rejecting  the  null  hypothesis 
that  there  is  no  damage.  The  only  situations  in  which  the  p-values  were  larger  were  for 
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1.5 
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0.5 
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01  23456789  10 

(a)  Randomly  perforated  damaged  geometry  for  (xd,  ho,  vd)  =  (3,  0.6,  0.3) 
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1.5 
1 

0.5 
0 

01  23456789  10 

(b)  Randomly  perforated  damaged  geometry  for  (xd,  hn,  vd)  =  (3,0.7,0.35) 

2 

1.5 

1 

0.5 
0 

01  23456789  10 

(c)  Randomly  perforated  damaged  geometry  for  (xq,  ho,  ho)  =  (3, 1.3, 0.65) 

2 

1.5 
1 

0.5 
0 

01  23456789  10 

(d)  Randomly  perforated  damaged  geometry  for  (xu,  hD,  vp)  =  (3, 1.4,  0.7) 

Figure  3:  The  example  perforated  damaged  geometry  for  (a)  (xd,  Hd^vd)  =  (3,  0.6,  0.3), 
(b)  for  (xD,hD,vD)  =  (3,0.7,0.35),  (c)  for  (xD,hD,vD)  =  (3,1.3,0.65),  (d)  for 
(xd,  ^DiVd)  =  (3, 1.4,  0.7),  respectively. 
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X  versus  area  using  left  nodes  and  xd  =  3 


X  versus  area  using  left  nodes  and  xd  —  5 
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nf3  P  versus  area  using  left  nodes  and  xd  =  3 
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p  versus  area  using  left  nodes  and  xd  =  5 
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Figure  4:  Results  using  the  left  nodes  for  observations: 
of  area  occupied  by  damage  for  (a)  xr,  =  3,  (b)  for  xr,  - 
versus  proportion  of  area  occupied  by  damage  for  xn  - 


(d) 

Test  statistic  X  versus  proportion 
=  5;  (c)  the  corresponding  p- values 
=  3,  and  (d)  for  xr,  =  5. 
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( ho,VD )  =  (0.6,  0.3)  and  ( Iid,vd )  =  (0.7,0.35)  so  we  will  focus  on  these  examples  in  the 
remainder  of  this  section. 

The  results  of  damage  detection  using  the  right  set  of  nodes  (see  Table  1)  are  depicted 
in  Figures  5(a)-(d).  We  see  similar  results  for  these  pixels  as  for  the  left  subset  of  pixels. 
In  the  example  where  xp  =  3,  both  ( hD,vn )  =  (0.6,  0.3)  and  ( hD,vn )  =  (0.7,0.35)  have 
high  p- values  when  calculated  using  Dtj  so  that  the  null  hypothesis  cannot  be  rejected  with 
even  50%  confidence  in  Figure  5(c).  Using  right  pixels  with  data  generated  with  xd  =  5, 
the  p-values  using  DtJ  are  low  for  both  ( hr>,VD )  =  (0.6,  0.3)  and  ( Iid,vd )  =  (0.7,0.35) 
with  rejection  of  the  null  hypothesis  with  90%  confidence  as  depicted  in  Figure  5  (d). 
The  result  of  this  test  using  DUnd  gives  a  low  p-value  for  the  example  damage  with 
(xd,  ho,  vd)  =  (5,  0.6,  0.3),  as  we  observe  in  Figure  5(d),  but  for  (5,0.7,0.35)  the  p-value 
is  over  0.25  and  we  may  only  reject  the  null  hypothesis  with  70%. 

We  turn  to  the  data  from  evenly  spaced  nodes  in  Figures  6  (c)  and  (d),  where  the  p- 
values  are  very  low  for  damage  over  3%.  For  the  smaller  damages  with  [hr,,  vr>)  =  (0.6,  0.3) 
and  (hD,vr)  =  (0.7,0.35),  the  p-values  calculated  using  D ^  are  high  so  that  one  cannot 
reject  the  null  hypothesis  with  very  high  confidence  as  can  be  seen  in  Figures  6  (c)  and 
(d).  The  p-values  using  D™nd  are  only  large  for  the  smallest  damage  in  Figure  6  (c)  and 
(d).  For  this  example  one  may  reject  the  null  hypothesis  with  95%  confidence  for  xd  =  3 
in  Figure  6(c)  and  85%  confidence  for  Xd  =  5  in  Figure  6(d). 

Finally,  in  Figures  7  (a)-(d)  we  depict  the  results  obtained  by  using  the  entire  set  of 
spatial  nodes  (all  seventeen  pixels).  All  examples,  except  the  smallest  example  damages 
[hr,  vd)  =  (0.6,  0.3)  for  both  xd  =  3  and  xd  =  5,  have  very  small  p-values  which  indicates 
that  the  null  hypothesis  using  all  spatial  nodes  can  be  rejected  for  damages  over  3%. 

In  summary,  our  simulations  suggest  that  for  the  damages  investigated  here,  the  meth¬ 
ods  we  propose  for  damage  detection  appear  useful  with  porous  materials. 
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X  versus  area  using  right  nodes  and  Xd  =  3 


X  versus  area  using  right  nodes  and  Xd  =  5 
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p  versus  area  using  right  nodes  and  x^,  =  3 


p  versus  area  using  right  nodes  and  xd  —  5 
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Figure  5:  Results  using  the  right  nodes 
of  area  occupied  by  damage  for  (a)  Xd 
versus  proportion  of  area  occupied  by 


for  observations:  Test  statistic  X  versus  proportion 
=  3,  (b)  for  xd  =  5;  (c)  the  corresponding  p-values 
damage  for  xd  =  3,  and  (d)  for  xd  =  5. 


13 


X  versus  area  using  even  nodes  and  Xd  =  3 


X  versus  area  using  even  nodes  and  Xd  =  5 
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p  versus  area  using  even  nodes  and  xd  =  3 


p  versus  area  using  even  nodes  and  xd  =  5 
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Figure  6:  Results  using  the  evenly  spaced  nodes  for  observations:  Test  statistic  X  versus 
proportion  of  area  occupied  by  damage  for  (a)  xd  =  3,  (b)  for  xd  =  5;  (c)  the  corre¬ 
sponding  p- values  versus  proportion  of  area  occupied  by  damage  for  xd  =  3,  and  (d)  for 
xd  =  5. 
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Figure  7:  Results  using  all  of  the  nodes 
of  area  occupied  by  damage  for  (a)  xr, 
versus  proportion  of  area  occupied  by 


for  observations:  Test  statistic  X  versus  proportion 
=  3,  (b)  for  xd  =  5;  (c)  the  corresponding  p- values 
damage  for  =  3,  and  (d)  for  =  5. 
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4  Damage  Characterization 

Though  the  detection  techniques  that  we  discussed  in  Section  3  can  be  used  to  detect 
damage,  we  would  also  like  to  characterize  the  extent  of  damage  once  its  presence  is 
suspected.  We  will  consider  the  ordinary  least  squares  (OLS)  estimation  procedure  in 
order  to  characterize  parameterized  damages.  Under  certain  conditions,  this  procedure 
can  provide  standard  error  estimates  for  the  damage  related  parameter  estimates.  The 
OLS  parameter  estimate  given  by 


q  =  argmin  J(q), 
q£Q 


(24) 


where  now  we  consider  the  parameter  set  q  =  (xd,  h,D,vn)  and  J  is  given  by  (16).  The 
OLS  error  variance  estimate  is  given  by 

£2  =  J(v) 

nm  —  p 

where  p  is  the  number  of  parameters  defining  the  dimension  of  the  admissible  parameter 
set  Q,  n  is  the  number  of  temporal  nodes  and  m  is  the  number  of  spatial  nodes  (pixels). 
The  covariance  matrix  estimate  [6,  13,  23]  is  then  given  by 


(25) 


£(?)  =  52  (x'(?)x(?))  *  . 

where  x(q)  is  the  matrix  of  sensitivities  with  entries 


_d_ 

9qk 


UiA  C) 


1 

C =<7 


(26) 


(27) 


for  i  =  1,2 , ,m,  j  =  1,2 , ,n,  and  k  =  1, 2, . .  .p.  The  details  of  the  calculations  of 
these  sensitivities  using  associated  sensitivity  partial  differential  equations  are  tedious  but 
straight  forward  and  will  not  be  given  here;  the  interested  reader  may  find  them  given  in 
Appendix  A. 2  of  [22],  The  asymptotic  standard  error  estimate  for  the  kth  parameter  is 
then  given  by 

SE (q)  «  y/xkk(q).  (28) 

We  will  use  a  very  large  number  of  nodes  so  for  all  our  parameter  sets,  the  95%  confidence 
interval  [6,  13]  is  given  by 


95%Confidence  Interval  =  (q  —  1.96  SE (q),  q  +  1.96  SE(g)) 


and  the  99%  confidence  interval  is  given  by 


99%Confidence  Interval  =  (q  —  2.58 SE (q),q  +  2.58 SE(g)). 

Unlike  the  situation  for  damage  detection,  there  are  significant  differences  in  the  results 
using  realizations  of  T)VJ  versus  realizations  of  DUnd.  All  three  variables  (x^,  hD,vr>) 
appear  to  be  well  estimated  using  the  realizations  of  D,;r  However,  there  are  many  cases 
where  using  realizations  of  DJ|nd  yield  values  that  are  significantly  different  than  the  true 
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Table  2:  The  estimated  geometries  for  (xo,  ho,vo)  =  (3, 1,0.5) 


Nodes 

n  x  m 

Data 

d2 

{xD,hD,vD) 

(SE(£D),SE(dD),SE(i)D)) 

Left 

1012 

D 

j-jrand 

0.028 

0.03 

(2.9,0.941,0.51) 

(3.84,0.653,0.59) 

(0.086,0.25,0.071) 

(0.073,0.18,0.099) 

Right 

1012 

D 

j-jrand 

0.027 

0.031 

(3.01,1.01,0.49) 

(4.09,0.626,0.58) 

(1.1,1.4,0.45) 

(0.084,0.2,0.12) 

Even 

1012 

D 

jprand 

0.028 

0.031 

(2.95,0.948,0.5) 

(3.88,0.564,0.64) 

(0.13,0.41,0.1) 
(0.083,  0.15,  0.12) 

All 

1564 

D 

jprand 

0.024 

0.027 

(2.93,0.934,0.51) 

(3.89,0.596,0.62) 

(0.08,0.26,0.07) 

(0.067,0.14,0.1) 

Table  3: 

The  estimated  geometries  for  (xo,  ho ,  ud)  =  (3,  2, 1) 

Nodes 

n  x  m 

Data 

d2 

(xD, hD,vD) 

(SE(xD),SE(dD),SE(uD)) 

Left 

1012 

D 

j^rand 

0.032 

0.048 

(3.01,1.99,1) 

(3.02,3.02,0.76) 

(0.023,0.059,0.015) 

(0.041,0.1,0.012) 

Right 

1012 

D 

j^rand 

0.03 

0.045 

(3.02,2.01,1) 

(3.6,1.71,0.93) 

(0.29,0.41,0.044) 

(0.27,0.48,0.072) 

Even 

1012 

D 

j-^rand 

0.03 

0.052 

(3.01,1.99,1) 

(3.02,2.87,0.78) 

(0.027,0.068,0.018) 

(0.044,0.11,0.017) 

All 

1564 

D 

j-jrand 

0.026 

0.06 

(3.01,1.98,1) 

(3.01,1.94,0.95) 

(0.021,0.054,0.015) 

(0.037,0.1,0.027) 

values.  In  some  of  these  cases,  the  standard  errors  are  sufficiently  large  so  that  a  99% 
confidence  interval  covers  these  values. 

In  Table  2,  which  presents  the  results  using  (xo,  ho,  Ud)  =  (3, 1,0.5),  the  estimate 
(xd,  hDl  vD)  for  the  data  D ^  is  much  closer  to  the  “true”  value  than  that  for  the  data 
DDnd.  For  this  example  using  a  realization  of  D)|nd,  for  each  set  of  pixels,  the  values 
( hr,,VD )  are  within  the  95%  confidence  intervals  though  the  estimated  values  are  not 
close  to  the  “true”  value.  The  estimates  of  Xo  for  each  set  of  pixels  are  not  only  far  from 
the  value  xo  =  3  but  the  standard  error  SE (xd)  is  sufficiently  small  so  that  for  each  pixel 
set  the  99%  confidence  interval  does  not  contain  the  value  xo  =  3. 

We  also  considered  the  realizations  and  D™nd  for  (xo,  ho,  Vo)  =  (3,2,1)  (see 
Table  3).  Again,  for  the  results  obtained  by  using  the  realization  D ij,  the  corresponding 
parameter  estimates  are  close  to  their  “true”  values  and  the  standard  errors  are  reasonable 
for  all  pixel  sets.  The  corresponding  realizations  of  D)|nd  with  all  pixel  sets,  except  the 
right  pixel  set,  yield  estimates  xo  that  are  very  close  to  the  value  xo  =  3.  For  every  set 
of  pixels,  the  value  xo  —  3  is  covered  by  the  95%  confidence  interval.  The  right  pixel 
set  and  the  full  pixel  set  (“All  17  Nodes”  in  Table  1)  have  values  of  ho  close  to  ho  =  2 
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Table  4:  The  estimated  geometries  for  (xo,  ho,vo)  =  (5, 1,0.5) 


Nodes 

n  x  m 

Data 

d2 

(xo,hD,VD) 

(SE(xD),SE(hD),SE(vD)) 

Left 

1012 

D 

j-jrand 

0.028 

0.032 

(5.08,1.04,0.49) 
(5.87, 1.12,  0.51) 

(0.1,0.47,0.086) 
(0.19, 0.65,  0.11) 

Right 

1012 

D 

j-jrand 

0.027 

0.031 

(4.92,0.944,0.51) 

(5.4,1.08,0.51) 

(0.088,0.26,0.071) 

(0.089,0.28,0.067) 

Even 

1012 

D 

jprand 

0.028 

0.033 

(4.95,0.934,0.52) 
(5.5, 1.12,  0.51) 

(0.12,0.37,0.094) 

(0.11,0.3,0.073) 

All 

1564 

D 

jprand 

0.024 

0.028 

(5.05,1.07,0.48) 

(5.64,1.39,0.45) 

(0.092,0.4,0.072) 

(0.12,0.41,0.052) 

Table  5: 

The  estimated  geometries  for  (xo,  ho,  vo )  =  (5,  2, 1) 

Nodes 

n  x  m 

Data 

d2 

(xD,  hD,  vo) 

(SE(xD),SE(dD),SE(dD)) 

Left 

1012 

D 

j^rand 

0.03 

0.038 

(4.99,2.01,1) 

(4.96,1.99,0.92) 

(0.037,0.098,0.02) 

(0.044,0.12,0.025) 

Right 

1012 

D 

j^rand 

0.029 

0.035 

(5,2,1) 

(5.17,2.31,0.86) 

(0.027,0.078,0.018) 

(0.04,0.1,0.017) 

Even 

1012 

D 

j-^rand 

0.03 

0.041 

(5,2.02,1) 

(5.12,2.22,0.87) 

(0.027,0.066,0.017) 

(0.04,0.1,0.02) 

All 

1564 

D 

j-jrand 

0.026 

0.038 

(4.99,2.01,1) 

(5.1,2.14,0.89) 

(0.022,0.058,0.015) 

(0.031,0.083,0.017) 

which  is  covered  by  their  respective  95%  confidence  intervals.  The  estimates  ho  for  the 
left  and  evenly  spaced  pixel  set  are  both  over  estimates  and  are  far  from  the  value  ho  =  2 
when  this  distance  is  compared  to  the  estimated  standard  errors.  Neither  the  left  nor 
evenly  spaced  subset  of  pixels  provided  estimates  in  which  ho  —  2  is  covered  by  the  95% 
confidence  interval.  We  see  that  the  estimates  Vo  in  this  case  are  close  to  the  value  Vo  —  1 
for  the  full  and  right  pixel  sets  while  the  estimates  Vo  are  much  smaller  than  Vo  =  1  for 
the  left  and  evenly  spaced  pixel  sets. 

The  results  of  the  OLS  parameter  estimation  procedure  for  realizations  of  D™nd 
with  (xo,  ho,  Vo)  =  (5, 1,0.5)  in  Table  4  are  similar  to  the  results  in  Table  2  with 
(xo,  ho,  vo)  =  (3, 1,0.5).  The  estimates  of  ( ho,VD )  are  close  to  ( ho,vo )  =  (1,0.5) 
and  the  standard  errors  (SE(Ad),  SE(Dd))  are  sufficiently  large  that  ( ho,i’o )  =  (1,0.5)  is 
within  two  standard  errors  of  the  parameter  estimate  ( Iio,vd )  for  all  sets  of  pixels  that 
we  considered.  The  estimate  xo  for  the  realization  D)|nd  is  much  larger  than  xo  =  5 
(though  this  effect  is  not  as  severe  as  for  xo  =  3  in  Table  2)  and  the  standard  errors  are 
small  so  that  for  each  set  of  pixels  the  value  xo  =  5  is  not  covered  by  the  resepective  95% 
confidence  interval. 

Using  the  realization  D™nd  in  Table  5  with  damage  given  by  (xo,  ho,  vo)  =  (5,  2, 1)  for 
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Table  6:  The  estimated  geometries  for  (xd, /id,  up)  =  (3,  0.5,1) 


Nodes 

n  x  m 

Data 

d2 

(■ xD,hD,vD ) 

{SE(xD),SE(hD),SE(vD)) 

Left 

1012 

D 

j-jrand 

0.028 

0.03 

(3.05,0.512,0.97) 

(3.39,0.987,0.66) 

(0.037,0.015,0.047) 

(0.062,0.16,0.063) 

Right 

1012 

D 

j^rand 

0.028 

0.032 

(2.88,0.389,1.2) 
(2.96,  0.476, 1.1) 

(0.066,0.087,0.51) 

(0.12,0.032,0.38) 

Even 

1012 

D 

j-jrand 

0.028 

0.031 

(2.99,0.513,0.97) 

(3.35,0.962,0.67) 

(0.047,0.018,0.068) 

(0.082,0.22,0.084) 

All 

1564 

D 

j-jrand 

0.025 

0.028 

(3.05,0.492,1) 

(3.41,0.954,0.67) 

(0.031,0.014,0.039) 

(0.061,0.14,0.061) 

each  estimate,  for  each  set  of  pixels,  the  estimates  are  within  0.22  of  the  “true”  parameter 
value.  Though  these  estimates  are  close  to  the  parameter  value,  many  of  the  standard 
error  estimates  are  still  so  small  that  the  “true”  parameter  value  is  not  covered  by  the 
99%  confidence  interval.  For  each  set  of  pixels,  the  99%  confidence  interval  of  i’o  does 
not  cover  the  value  vp  =  1. 

We  also  considered  damages  defined  through  ellipses  that  have  a  vertical  semi-axis 
which  is  larger  than  the  horizontal  semi-axis.  In  these  examples,  this  inequality  was 
not  maintained  in  the  parameter  estimates  which  were  computed  using  realizations  of 
D™nd.  In  Table  6,  we  present  the  results  for  (iD,  hD,  vd)  =  (3,  0.5,1).  The  estimates  of 
(xD,hD,V]j)  for  the  realization  D jj  are  farther  from  the  “true”  value  than  in  the  previous 
examples.  This  implies  that  the  geometry  of  these  ellipses  poses  a  problem  in  the  param¬ 
eter  estimation  problem.  This  effect  is  amplified  for  the  realizations  DDnd.  For  the  pixel 
sets  “Left,”  “Even,”  and  “All”  the  parameter  estimation  procedure  produced  an  ellipse 
estimate  with  horizontal  semi-axes  longer  than  vertical  semi-axes.  Also,  in  these  exam¬ 
ples  the  estimate  of  the  center  xd  is  larger  than  the  value  xd  =  3.  The  pixel  set  “Right” 
produced  an  estimate  which  was  very  close  to  the  “true”  value  with  each  95%  confidence 
interval  containing  that  “true”  parameter  value.  The  estimates  using  the  realization  of  Do¬ 
using  the  “Right”  pixel  set  is  the  farthest  from  the  true  damage  ( xd ,  hp,  vd )  =  (3,  0.5, 1) 
while  the  parameter  estimate  using  the  “Right”  with  realizations  of  D)“ld  is  the  closest. 
This  effect  occurs  because  for  data  generated  without  pores,  the  “Right”  pixels  are  the 
farthest  from  the  damage  centered  at  xp  —  3,  so  there  is  relatively  little  information  in 
those  pixel  sets  compared  to  the  other  pixel  sets.  For  data  generated  using  realizations 
of  D)|nd,  the  smaller  the  distance  from  the  bottom  of  the  damage  to  the  source  boundary 
(uq),  the  larger  the  error  associated  with  (uv^nd)ij  so  the  “Right”  pixel  set  contains  the 
least  amount  of  error  clue  to  the  approximation  of  {ur^nd)ij  with  Uij.  We  see  very  similar 
results  for  the  damage  (xo,  hp,  vp)  =  (3,  0.5, 1.5)  in  Table  7. 

When  a  damage  with  vd  >  h d  is  placed  in  the  horizontal  center  of  the  10  by  2 
rectangle,  all  of  the  pixel  sets  contain  pixels  directly  below  the  center  of  the  damage. 
Realizations  of  D;/  and  l);)"'1  both  yield  relatively  small  standard  errors  compared  to 
the  difference  between  the  parameter  estimate  and  the  true  parameter  set.  It  appears 
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Table  7:  The  estimated  geometries  for  (i0,  hp,  vp)  =  (3,  0.5, 1.5) 


Nodes 

n  x  m 

Data 

d2 

(xd,  h.D,  VD ) 

{SE(xD),SE{hD),SE{vD)) 

Left 

1012 

D 

j-jrand 

0.029 

0.043 

(2.97,0.5,1.5) 

(3.14,1.1,0.9) 

(3.9  x  10~3,  6.8  x  10“3,  0.015) 
(0.039,0.099,0.057) 

Right 

1012 

D 

j-jrand 

0.028 

0.037 

(2.97,0.469,1.6) 

(2.59,0.546,1.8) 

(0.015,0.012,0.06) 

(0.051,0.045,0.22) 

Even 

1012 

D 

j-jrand 

■ 

(2.95,0.5,1.5) 
(3.19, 1.1,0.91) 

(5.5  x  10~3,  8.3  x  10“3,  0.017) 
(0.06,0.14,0.069) 

All 

1564 

D 

j-jrand 

0.025 

0.035 

(2.97,0.506,1.5) 

(3.18,1.14,0.88) 

(3.6  x  10-3, 6.4  x  10~3,  0.015) 
(0.039,0.1,0.051) 

Table  8:  The  estimated  geometries  for  (xD,hD,vo)  =  (5,  0.5,1) 


Nodes 

n  x  m 

Data 

d2 

(xD,hD,i)D) 

(SE(xD),SE(M,SE(t)D)) 

Left 

1012 

D 

j-jrand 

(4.98,0.447,1.1) 
(6.18,  0.633, 1.1) 

(0.021,0.016,0.036) 

(0.061,0.061,0.07) 

Right 

1012 

D 

j^rand 

0.028 

0.036 

(4.98,0.439,1.1) 

(4.87,0.729,0.86) 

(0.023,0.016,0.039) 

(0.051,0.065,0.065) 

Even 

1012 

D 

j-jrand 

(4.99,0.47,1) 
(5.31, 1.2, 0.62) 

(0.031,0.016,0.051) 

(0.086,0.21,0.061) 

All 

1564 

D 

j-jrand 

(4.98,0.483,1) 
(5.29,  3.18,  0.41) 

(0.027,0.012,0.042) 

(0.13,0.44,0.016) 

20 


Table  9:  The  estimated  geometries  for  (xd,  Ad,Ld)  =  (5,  0.5, 1.5) 


Nodes 

n  x  m 

Data 

d2 

(xD,hD,vD) 

{SE(xD),SE{hD),SE(vD)) 

Left 

1012 

D 

j-jrand 

0.028 

0.054 

(5.03,0.482,1.6) 

(5.92,0.641,1.7) 

(2.9  x  10~3,  6.4  x  10"3, 0.012) 
(0.014,6.9  x  10“3,  0.039) 

Right 

1012 

D 

j-jrand 

0.029 

0.053 

(4.97,0.487,1.5) 

(4.75,0.834,1.2) 

(4.7  x  10~3,  6.4  x  10“3,  0.014) 
(0.035,0.032,0.047) 

Even 

1012 

D 

j-jrand 

0.028 

0.045 

(5,0.491,1.5) 

(5.15,1.28,0.84) 

(4  x  10~3,  7.3  x  10”3, 0.017) 
(0.071,0.14,0.053) 

All 

1564 

D 

j-jrand 

0.025 

0.039 

(5,0.491,1.5) 

(5.17,1.42,0.78) 

(3.2  x  10-3, 5.9  x  10~3,  0.013) 
(0.044,0.12,0.039) 

that  in  this  case,  the  damage  is  not  as  well  characterized.  In  Table  9,  with  (xd,  ho,  vd )  = 
(5,  0.5, 1),  the  parameter  estimates  for  the  realization  T)VJ  are  closer  to  (5, 0.5,1)  than  those 
of  the  realization  D™nd.  In  many  of  these  cases,  the  95%  confidence  interval  does  not 
contain  the  true  parameter  value  for  this  realization  of  D Using  the  realization  D™nd 
in  this  example,  the  parameter  estimates  with  pixel  sets  “Left”  and  “Right”  maintain 
I’d  >  ho  while  parameter  estimates  with  pixel  sets  “Even”  and  “All”  invert  this  inequality 
resulting  in  the  incorrect  characterization  vd  <  ho-  The  estimate  of  xd  =  6.18  using  the 
“Left”  pixel  set  for  the  realization  DDnd  is  very  far  from  xd  —  5  while  the  estimates 
of  xd  for  pixel  sets  “Right,”  “Even”  and  “All”  are  within  0.31  of  xd  =  5.  Though 
these  estimates  are  different  from  the  “true”  parameter  value,  the  standard  errors  are 
relatively  small,  resulting  in  the  “true”  parameter  value  being  many  standard  errors  from 
the  parameter  estimate.  Again,  for  the  larger  damage  (xd,1id,vd)  =  (5,  0.5, 1.5),  we  see 
similar  results  in  Table  9. 

In  summary,  our  approximate  homogenization  model  appears  to  work  very  well  to 
detect  damage  while  to  characterize  damage  the  homogenization  approximation  works  well 
only  in  some  cases  but  also  is  somewhat  less  satisfactory  in  other  cases.  This  is  dependent 
on  shape  of  the  damage  (it  appears  to  characterize  damage  with  Ad  >  Ld  better  than 
damage  with  ho  <  vd)  and  the  choice  of  pixel  set.  One  possible  reason  behind  some 
of  these  inaccuracies  in  characterizing  damage  may  be  attributed  to  the  large  residuals 
obtained  in  many  cases.  We  observe  that  large  residuals  produce  larger  values  of  the 
estimated  variances  a2  (see  (25),  (26),  and  (28))  which  in  almost  all  table  values  are  an 
order  magnitude  too  large  (recall  cr0  =  0.015  in  all  data  and  D™nd  generated).  This 
leads  to  large  SE  and  hence  larger  confidence  intervals.  This  could,  if  these  are  given 
too  much  credence,  lead  to  less  than  accurate  assessments  of  one’s  ability  to  characterize 
the  damage  with  these  methods.  Finally,  for  the  given  geometry  investigated  here  using 
thermal  probes,  it  is  not  surprising  that  characterization  for  damage  with  Vo  >  ho  is  not 
as  accurate  as  that  for  examples  with  ho  >  vd-  We  note  that  the  latter  damages  are 
expected  to  be  more  prevalent  in  composite  materials  of  the  type  considered  here  where 
the  damages  tend  to  develop  from  enlarged  “horizontal”  pores  resulting  from  the  layered 
manufacturing  process  [3,  4,  5]. 
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5  Conclusions 


Based  on  previous  results  in  [3,  4,  5],  we  considered  using  parameter  estimation  proce¬ 
dures  in  homogenized  models  to  detect  and  characterize  damage  in  porous  materials.  We 
found  that  these  techniques  could  be  used  to  detect  damage  very  well  for  the  considered 
damages  in  porous  materials  which  were  above  3%  of  the  domain’s  area.  In  these  ex¬ 
amples,  there  was  no  clear  difference  between  using  data  generated  using  the  solution  of 
the  heat  equation  on  a  random  perforated  domain,  and  data  generated  using  the  homog¬ 
enization  approximation.  On  the  other  hand,  in  characterization  of  elliptical  damage, 
there  does  appear  to  be  a  difference  in  using  data  generated  using  the  solution  of  the  heat 
equation  on  a  perforated  domain  rather  than  using  the  homogenization  approximation 
in  estimating  the  center  and  the  size  of  the  elliptical  damage.  As  one  would  expect,  the 
estimates  of  the  center  and  semi-axes  lengths  using  data  generated  using  the  homoge¬ 
nization  approximation  were  always  closer  to  the  true  value  than  estimates  of  the  center 
and  semi-axes  lengths  using  data  generated  using  the  solution  of  the  heat  equation  on 
a  perforated  domain.  The  accuracy  of  the  estimates  of  the  center  coordinate,  horizon¬ 
tal  semi-axis  and  the  vertical  semi-axis  (xo,  ho,  vjf)  varied  over  the  examples  using  data 
generated  the  solution  of  the  heat  equation  on  the  randomly  perforated  domain.  The 
examples  with  hD  <  Vo  did  not  preserve  this  inequality  in  the  estimates  of  ho  and  Vo  in 
most  cases  while  in  the  examples  with  ho  >  Vfj  the  inequality  was  preserved.  When  we 
estimated  (xo,  ho,i’o)  using  data  generated  using  the  solution  of  the  heat  equation  on 
the  perforated  domain,  we  observed  significant  differences  between  using  different  pixel 
sets.  It  appeared  that  some  of  the  pixel  sets  were  more  sensitive  to  the  error  associated 
with  using  the  homogenization  approximation  in  the  inverse  problem.  It  is  worthwhile  to 
note  that  in  most  cases  the  estimated  area  of  the  damage  was  larger  in  the  estimate  than 
the  area  of  the  “true”  damage.  This  suggests  that  this  estimation  procedure  may  not 
estimate  the  center  nor  the  size  correctly,  often  yielding  an  estimate  that  is  larger.  These 
methods  could  be  used  as  a  first  estimate  of  damage  which  could  be  further  resolved  using 
a  subsequent  nondestructive  evaluation  technique. 
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